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TROOUCTION 


The  accurate  prediction  or  experimental  determination  of  projectile 
aerodynamics  Is  of  significant  importance  to  the  shell  designer  and  ballis- 
ticlan.  The  shell  designer  requires  accurate  aerodynamic  data  for  the 
overall  development  of  new  shell.  The  ballisticlan  Is  concerned  with  the 
development  of  aiming  data  and  therefore  relies  heavily  on  accurate  aero¬ 
dynamic  data.  Experimental  costs  have  sky  rocketed  In  recent  years  and 
have  contributed  significantly  to  overall  system  development  costs.  Compu¬ 
tational  techniques  are  beginning  to  show  promise  as  a  means  to  alleviate 
or  at  least  temper  these  rising  development  costs  by  providing  relatively 
low  cost  computer  analysis  of  new  designs.  As  computer  technology^ 
Increases  and  machines  become  faster  with  larger  memory,  the  use  of  compu¬ 
tational  methods  In  design  becomes  more  of  a  reality.  _ _ 


The  means  to  compute  projectile  aerodynamics  for  all  Mach  number 
regimes  covered  by  a  given  projectile  In  its  flight  history  has  been  an 
area  of  research  actively  pursued  by  the  Aerodynamics  Research  Branch  of 
the  Ballistic  Research  Laboratory.  Early  work  had  focused  on  the  super¬ 
sonic  flight  regime  and.  In  particular  on  the  accurate  prediction  of  the 
Magnus  force.  The  Magnus  force,  which  Is  very  small  In  magnitude  (on  the 
order  of  1/10  the  normal  force).  Is  a  critical  parameter  In  determining  the 
dynamic  stability  of  shell.  The  Magnus  force  Is  generated  by  a  spin 
Induced  distortion  of  the  boundary  layer;  therefore,  correct  modeling  of 
the  vlscous/lnvlscld  Interatlon  Is  critical  for  accurate  computations.  The 
work  of  Sturek,  et.  al.1  has  shown  that  accurate  results  In  the  supersonic 
regime  can  be  obtained  for  ogive- cylinder  projectile  shapes.  This  techni¬ 
que  Involved  separate  computations  of  the  turbulent,  viscous  boundary  layer 
and  the  outer  Invlscld  flow  field.  As  the  projectile  shapes  were  general¬ 
ized  to  Include  boattalls,  more  sophisticated  computational  techniques  had 
to  be  employed.  These  new  methods,  which  solved  the  thin-layer  Navler 
Stokes  equations,  were  successfully  applied  to  oglve-cyllnder-boattall 
shapes  by  Sturek  and  Schlff2'3.  The  solution  of  the  Navler  Stokes 
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equations  allows  for  the  simultaneous  computation  of  the  viscous/ In vise Id 
flow  field  and  thus  provided  the  basis  for  good  Magnus  force  prediction. 

A  region  of  critical  aerodynamic  behavior  occurs  in  the  high  transonic 
regime,  0.90  <  M  <  1.15,  where  aerodynamic  coefficients  have  been  found  to 
Increase  by  as  much  as  100%.  This  flight  velocity  regime  Is  both  experi¬ 
mentally  and  computationally  difficult.  Thus,  only  a  small  amount  of 
experimental  data  are  available  for  design  studies  and  only  limited  compu¬ 
tational  studies  have  been  made.  An  Initial  attempt  to  develop  a  computa¬ 
tional  capability  suitable  for  Nignus  prediction  at  transonic  velocity  was 
made  by  Nletublcz,  et.  al4  whereby  the  thln-layef  Navler  Stokes  computa¬ 
tional  technique  was  applied  to  standard  and  ho/low  projectile  shapes  at 
zero  degrees  angle  of  attack.  Computational  rekults  were  also  obtained  for 
non-spinning  projectiles  at  angle  of  attack s;  hUrfsver  limitations  of 
computer  resources  (COC  7600)  became  apparent.  "^In  this  paper,  recent 
results  are  presented  for  Magnus  force  computations  using  the  Cray  IS 
computer.  Comparisons  are  made  with  the  earlier  COC  computer  results  and 
are  further  compared  to  some  limited  experimental  data.  A  discussion  of 
the  numerical  technique  Is  Included.  These  results  represent  the  first 
computations  of  the  Magnus  effect  on  projectile  shapes  In  the  transonic 
flight  regime. 

II.  GOVERNING  EQUATIONS  AND  COMPUTATIONAL  TECHNIQUE 

The  set  of  equations  which  govern  fluid  motion  are  the  Navler-Stokes 
equations.  This  set  of  highly  non-linear  partial  differential  equations 
have  proven  extremely  difficult  to  solve  In  their  entirety.  In  most 
Instances  approximations  had  to  be  made  before  a  solution  could  be  attemp¬ 
ted.  For  example,  the  well  known  boundary  layer  equations  are  derived  by 
applying  approximations  to  the  Navler  Stokes  equations.  The  present  solu 
tlon  technique  also  makes  use  of  an  approximation.  The  Navler-Stokes  equa¬ 
tions  solved  here  make  use  of  the  thin-layer  approximation.  That  Is,  the 
viscous  terms  are  neglected  In  both  the  longitudinal  and  circumferential 
directions.  The  viscous  terms  are  retained,  however.  In  a  direction  nearly 
normal  to  the  surface  where  large  flow  field  gradients  exist.  This  formu¬ 
lation  retains  the  momentum  equations  In  all  three  coordinate  directions. 
The  retention  of  the  three  momentum  equations  allows  for  the  computation  of 
separated  flow  and  thus  differs  significantly  from  boundary  layer 
assumptions. 


The  equations  solved  here  are  written  In  a  generalized  coordinate 
system.  This  allows  a  wide  variety  of  body  shapes  to  be  computed  using  the 
same  basic  numerical  technique.  The  notation  for  the  physical  and  trans¬ 
formed  coordinate  systems  are  shown  In  Figure  1.  The  three  dimensional, 

transformed,  thin-layer  Navler-Stokes  equations,  written  In  non-  _ 

dimensional,  strong  conservation  law  form  are6 
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The  general  coordinate  transformations  are 
defined  as 
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A 

The  viscous  matrix,  S,  Is  written  as 


l-l 


0 

w(  ( u/3)  ( Cxu?+  C/5+c2w5)  «x 

uUj+C2*?2)  V;+  Cli/3)  («xV5y  Wc^y 
u(cj*-cj+cj)w5+(u/3)(cxuc+cyvc+c2wc)c2 
{ ( C  ^  C  J)  Cu/2(  u2*  v2*  w2)  5+ rf»r‘ 1  ( Y-1 )  * 1  ( «  2) 

♦ ( m/3) ( Cxu+  Cy v+c2w) ( cxuc+  Cyv?+  c2wc) } 


The  velocities 


U  ■  5t  +  5xu  +  £yv  + 

V  «  nt  +  nxu  +  nyv  +  n2w  (2) 

M  "  ct  +  5xu  +  cyV  +  5z" 

represent  the  contravarlent  velocities.  The  non-dimensional  velocities  U, 

V,  an  W  are  those  components  In  the  direction  of  the  transformed  coordi¬ 
nates  E,  n,  and  5,  respectively.  The  Cartesian  velocity  components  u,  v,  w 
together  with  the  density  p  and  total  energy  per  unit  volume  e  are  retained 
as  the  dependent  variables.  The  local  pressure,  p.  Is  determined  using  the 
relation 


p  •  (y  -  l)(e  -  .5p(u2  ♦  v2  +  w2)). 


The  velocities  are  non-dlmenslonallzed  by  the  free  stream  speed  of  sound 
a,,  the  density  by  p„,  and  the  total  energy  by  P.a2.  The  additional 

parameters  appearing  In  equation  1  are:  (a)  coefficient  of  thermal  con¬ 
ductivity,  k;  (b)  dynamic  viscosity,  u;  (c)  Reynolds  number  based  on  body 
diameter.  Re;  (d)  Prandtl  number,  Pr,  (e)  X  which,  based  on  Stokes 
hypothesis.  Is  -2/3  ». 
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As  mentioned  earlier,  these  equations  are  written  In  transformed  coord¬ 
inates;  therefore,  the  various  body  shapes  are  Introduced  through  determi¬ 
nation  of  the  metric  terms  5X,  nx,  cz,  etc..  These  terms  are  formed  by  a 

combination  of  the  derivative  terms  x?,  y?,  z?,  etc.,  and,  together  with 

the  transformation  Jacobian,  allow  for  variable  body  geometries.  Thus,  one 
of  the  first  steps  In  performing  a  computation  fs  the  generation  of  a 
computational  grid  which  provides  the  x,  y,  z  points  for  the  metric  deter¬ 
mination.  These  points  are  determined  prior  to  the  computations  and  are 
not  changed  with  time.  Examples  of  the  computational  grid  used  In  this 
study  are  shown  In  Figures  2  and  3.  A  two  dimensional  slice  of  the  overall 


grid  Is  shown  In  Figure  2.  The  upstream,  downstream  and  outer  flow  field 
computational  boundaries  extended  approximately  18  body  diameters  from  the 
body  surface.  At  this  distance  the  flow  field  should  be  uniform  and  the 
Imposed  boundary  conditions  are  considered  valid.  Figure  3  shows  an 
expanded  view  near  the  body.  The  clustering  of  grid  points  near  the  body 
surface  Is  required  In  order  to  resolve  the  viscous  components  of  the  flow 
field  near  the  body  surface.  Due  to  the  lack  of  sufficient  computer 
storage,  judicious  use  must  be  made  of  the  limited  grid  points  available. 

In  regions  where  the  viscous  effects  are  not  predominant  and  the  flow  field 
changes  slowly,  the  grid  points  are  sparse.  Additional  grid  clustering  Is 
used  In  the  longitudinal  direction  where  flow  field  gradients  are  expected. 
The  two  dimensional  grid  shown  In  Figure  2  was  rotated  about  the  body  axis 
In  10  degree  Increments  In  order  to  obtain  the  three  dimensional  grid 
required  for  computations  at  angle  of  attack. 
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As  mentioned  In  the  Introduction,  the  Magnus  effect  Is  produced  by  a 
spin  Induced  distortion  of  the  boundary  layer.  The  computation  must  there¬ 
fore  be  fully  three  dimensional  since  no  plane  of  symmetry  exists.  The 
boundary  conditions  used  for  the  computations  are: 

(1)  Inner  boundary,  body  surface 

U  •  W  *  0 

V  »  m  ■  non-dimensional  angular  velocity 

p  -  first  order  extrapolation 

p  -  calculated  using  the  three  transformed  momentum 
equations 

(2)  outer  boundary 

constant  free  stream  values  are  used  for  all  variables 

(3)  downstream  boundary 

M  <  1  pressure  Is  fixed  at  pM  and  all  other  variables 
are  extrapolated 

M  >  1  first  order  extrapolation  on  all  variables 

The  numerical  scheme  used  for  the  solution  of  equation  1  Is  a  fully 
Implicit,  approximately  factored,  finite  difference  algorithm  In  delta  form 
as  analysed  by  Beam  and  Warming7.  This  scheme  can  be  first  or  second  order 
accurate  In  time  and  second  or  fourth  order  accurate  In  the  three  spatial 
directions.  The  solution  of  the  three  dimensional  equations  Is  Implemented 
by  an  approximate  factorization  which  allows  the  system  of  equations  to  be 
solved  In  three  coupled  one  dimensional  steps.  This  procedure  has  been 
utilized  In  previous  applications2’5  with  a  high  degree  of  success.  Addi¬ 
tional  details  of  the  numerical  method,  computational  algorithm  and 
boundary  conditions  can  be  found  In  Reference  6. 

The  turbulence  model  employed  Is  an  algebraic  eddy  viscosity  model  as 
developed  by  Baldwin  and  Lomax5.  This  same  model  has  yielded  excellent 
results  for  Magnus  effects  at  supersonic  velocities3. 

As  Indicated  In  equation  1,  this  solution  technique  Involves  solving 
the  time  dependent  Navler-Stokes  equations.  The  procedure  Is  started  by 
assuming  a  uniform,  free  stream  solution  for  all  grid  points  In  the  compu¬ 
tational  domain.  The  calculation  then  marches  In  time  until  a  steady  state 
solution  Is  obtained.  The  Implicit  technique  used  here  allows  for  large 
time  steps  to  be  taken  which  helps  to  reduce  the  computation  time. 
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III.  MOOEL  GEOMETRY  AND  EXPERIMENT 

One  means  c*  establishing  the  computational  accuracy  of  a  numerical 
scheme  Is  through  comparisons  with  available  experimental  data.  The  model 
used  for  the  experiment  and  computational  study  presented  here  Is  an 

Idealization  of  a  realistic 
artillery  projectile 
geometry.  The  experimental 
model  shown  In  Figure  4 
consists  of  a  3  caliber  (1 
caliber  *  maximum  body 
diameter)  sharp  secant-ogive 
nose,  a  two  caliber  cylin¬ 
drical  mid-section,  and  a  1 
caliber  7°  conical  afterbody 
or  boattall.  A  similar 
model  was  used  for  the 
computational  studies  with 
the  only  difference  being  a 
5%  rounding  of  the  nose  tip. 
The  nose  tip  rounding  was  done  for  computational  efficiency  and  Is  consi¬ 
dered  to  have  little  Impact  on  the  final  Integrated  forces. 

The  surface  pressure  experimental  data9  used  for  comparison  In  this 
paper  was  obtained  In  the  NASA  Langley  8  Ft.  Pressure  Tunnel.  The  test 
conditions  of  1  atmosphere  supply  pressure  and  320°K  supply  temperature 
resulted  In  a  Reynolds  number  of  4.5  x  10s  based  on  model  length.  The 
model  was  Instrumented  with  pressure  ports  at  15  longitudinal  stations. 
Pressure  data  and  aerodynamic  force  data  were  obtained  at  Mach  numbers  of 
0.91  to  1.2  and  angles  of  attack  from  0  to  10  degrees. 

Additional  tests  were  conducted  at  the  Naval  Surface  Weapons  Center  for 
similar  tunnel  conditions.  This  test  utilized  the  same  model  which  now 
Included  the  capability  for  spin.  The  model  was  spun  to  500  rev/ sec;  and 
Magnus  force  and  moment  measurements  were  obtained  while  the  model  coasted 
down  to  zero  spin.  Aerodynamic  force  measurements  were  obtained  at  Mach  * 
0.91  for  a  «  O^  to  10°. 


IV.  COMPUTER  RESOURCES 

The  Initial  computations  were  carried  out  on  the  BRL  CDC  7600  computer. 
This  machine  has  a  maximum  large  core  storage  capability  of  approximately 
380,000  useable  words.  This  limited  the  maximum  computational  grid  to  be 
21,600  points  (60  longitudinal,  20  normal  and  18  circumferential).  As 
previously  mentioned  the  solution  Is  marched  In  time  until  a  steady  state 
condition  is  obtained.  A  typical  converged  solution  required  approximately 
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1000  time  steps.  The  COC  computation  ran  at  a  speed  of  12.96  x  10~4 
sec/time  step/point.  This  resulted  In  7.78  hour  computation  time  for  one 
converged  solution  at  one  set  of  conditions.  As  will  be  shown  later  this 
computational  grid  did  not  provide  adequate  grid  resolution. 

In  order  to  obtain  Improved  grid  resolution  and  faster  computational 
speed,  the  computer  code  was  placed  on  the  Cray  IS  computer  at  Klrtland  Air 
Force  Base.  Operationally,  the  connection  to  the  Cray  was  made  via  the 
ARPANET  from  Aberdeen  Proving  Ground.  The  ARPANET  Is  a  Department  of 
Defense  digital  switching  network  which  allows  terminals  and  geographically 
separated  computers  to  communicate.  The  Cray  IS  computer  Is  a  vector 
processor  and  has  a  demonstrated  speed  advantage  over  the  CDC  7600  of  10  In 
the  vector  mode.  The  present  computations  have  not  yet  taken  full  advan¬ 
tage  of  the  vector  capability,  however.  The  same  computation  described 
earlier  ran  at  a  speed  of  5.2  x  10“4  sec/time  step/point  on  the  Cray  which 
resulted  In  a  2.5x  Increase  In  speed  for  the  same  grid.  An  additional 
advantage  in  using  the  Cray  Is  the  Increased  storage  capability.  The 
present  configuration  allows  for  1  million  words  of  storage  with  an  expan¬ 
sion  capability  to  4  million  words.  The  computational  grid  was  therefore 
Increased  to  51,840  grid  points  (60  longitudinal,  24  normal  and  36  circum¬ 
ferential).  It  was  felt  that  Increased  grid  resolution  in  the  circumfer¬ 
ential  direction  was  required  for  determination  of  the  Magnus  force  compon¬ 
ents.  With  the  Increased  grid  and  same  convergence  crlterial  used  earlier, 
the  final  computational  time  was  7.49  hours  on  the  Cray.  With  full  vector- 
Izatlon,  It  Is  anticipated  that  this  time  can  be  reduced  to  1.5  hours. 

With  the  advent  of  even  faster  machines  (Cray  2,  with  32  million  words  of 
storage  and  speeds  6-12  times  faster  than  the  Cray  1)  the  eventual  computa¬ 
tional  time  can  potentially  be  reduced  to  15  minutes. 


V.  RESULTS 

The  results  to  be  presented  are  In  the  form  of  pressure  coefficients, 
velocity  profiles  and  aerodynamic  coefficients.  The  computed  results  are 
compared  with  the  experimental  data  for  M  *  0.91,  a  »  2.0°,  and  Re  * 

4.5  x  io6.  The  Magnus  measurements  were  compared  at  PO/U,.  ■  0.39  where  P  * 
angular  velocity,  D  ■  maximum  body  diameter  and  U,,  *  free  stream  velocity. 
Surface  pressure  comparisons  are  made  using  the  Langley  data,  while  the 
aerodynamic  coefficient  comparisons  are  made  with  the  NSWC  data10. 

a.  Surface  Pressure  Coefficient 

The  surface  pressure  coefficient,  Cp  ■  (p  -  pj/.5pju£,  on  the  leeward 

ray  is  shown  as  a  function  of  longitudinal  position  In  Figure  5.  The 
experimental  data,  and  computational  results  from  both  the  CDC  7600  and 
Cray  IS  computers  are  shown.  The  CDC  results  (dashed  line)  show  marginal 
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agreement  with  the  experimental  data  (circles)  over  the  entire  projectile 

surface.  The  expansions  and 
recompressions  which  occur  near 
the  ogive- cylinder  and  cyllnder- 
boattall  junction  have  not  been 
adequately  captured  by  the  compu¬ 
tation.  The  solid  line  in  Figure 
5  Is  the  results  from  the  compu¬ 
tation  on  the  Cray  IS  computer 
where  the  grid  has  been  expanded 
to  36  points  circumferentially, 

24  points  normal  to  the  surface 
and  the  60  longitudinal  points 
have  been  redistributed.  The 
agreement  with  the  experimental 
data  has  Improved  significantly. 
Some  discrepancy  Is  still  appar¬ 
ent  however,  on  the  cylinder  and 
on  the  boattail.  As  can  be  seen 
from  Figure  3  the  grid  points  In 
the  vicinity  of  the  boattail  have  been  severely  stretched  and  Is  the  appar¬ 
ent  cause  of  the  discrepancy.  All  remaining  results  to  be  presented  were 
obtained  from  the  Cray  IS.  A  comparison  of  the  windward  and  leeward  pres¬ 
sure  distribution  Is  shown  In  Figure  6.  The  ogive  experiences  high  pres¬ 
sure  along  the  windward  side 
whereas  the  high  pressure  for  the 
boattail  Is  on  the  leeward  side. 
On  the  ogive,  the  high  pressure 
on  the  windward  ray  causes  an  up¬ 
ward  force.  On  the  boattail,  the 
high  pressure  on  the  leeward  ray 
causes  a  negative  or  downward 
force.  This  condition  forms  a 
couple  about  the  center  of 
gravity  and  contributes  to  the 
critical  aerodynamic  behaviour 
which  occurs  at  transonic 
velocities. 


Surface  Pressure  Coefficient,  Exoerinent  and 
Canputatlms.  f*  «  0.91,  a  •  2.0* 


Figure  9.  Ulnduerd  and  Leeward  Surface  Pressure 
Coefficient,  X  ■  0.91,  a  •  2.0* 
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Comparisons  between  computation  and  experiment  for  the  circumferential 
surface  pressure  distribution  are  shown  In  Figures  7  and  8  for  X/D  stations 
1.56  and  5.19  respectively.  As  seen  In  these  figures*  the  computation 
predicts  the  correct  trend  of  the  data  while  the  actual  magnitude  shows  a 

deviation.  Upon  close  evaluation 
of  the  circumferential  pressure 
distribution,  a  slight  asymmetry 
Is  detected  about  the  +  «  180 
plane.  This  condition  contri¬ 
butes  to  the  expected  side  force. 
The  grid  used  for  these  calcula¬ 
tions  Is  by  no  means  optimum  and 
additional  computational  experi¬ 
mentation  Is  required.  The 
results  do,  however.  Indicate  the 
potential  for  obtaining  satisfac¬ 
tory  aerodynamic  coefficients. 
This  optimism  Is  based  on  results 
which  have  been  obtained5 
previously  for  an  axl symmetric 
cqpputatlon.  Figure  9  shows  a 
comparison  between  experimental 
data  and  a  computation  using  a  grid  consisting  of  80  longitudinal  points 
and  40  points  In  the  near  normal  direction.  The  agreement  Is  excellent  and 
Indicates  the  quality  of  results  which  should  be  possible  for  the  three 
dimensional  cases  given  adequate  grid  resolution. 


Figure  7.  Clrcunferentlel  Pressure  Distribution, 
M  -  0,91,  a  «  2.0*.  X/0  *  1.S6 


figure  I.  Ctrcueferentlel  Pressure  Distribution,  F,9«re  Surface  Pressure  Coefficient,  Ixperiaent 

H  •  0.91,  o  •  2.0*,  X/0  ■  S.lt  Axl syene trie  Caneutitlon,  *  •  0.9*.  a  •  0.0* 
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A  primary  purpose  of  this  research  effort  Is  the  development  of  a  capa¬ 
bility  for  the  prediction  of  aerodynamic  coefficients  and.  In  particular, 
the  capability  to  compute  the  Magnus  effect.  As  noted,  the  Magnus  effect 
Is  a  viscous  phenonena  associated  with  the  spinning  projectile.  Therefore, 
In  order  for  a  computational  technique  to  predict  this  effect.  It  must 
adequately  compute  the  longitudinal  and  circumferential  wall  shear  stress 
for  the  spinning  projectile  at  angle  of  attack.  The  experimental  determin¬ 
ation  of  the  u,  v,  and  w  velocity  distribution  is  especially  difficult  at 
transonic  velocities.  Although  no  experimental  data  are  available  for 
comparison,  the  computed  circumferential  velocity  distributions  are  shown 
In  Figures  10  and  11  for  X/D  *  4.22  and  X/D  «  5.50,  respectively.  A  signi¬ 
ficant  asymmetry  can  be  seen  In  the  velocity  distributions  at  t  »  90°  and 


Figure  1*1.  Circumferential  Velocity  Profile!  for 
9*0*.  VO* ,  180* ,  270*  at  *  •  0.9i, 
a  *  2.0”,  X/D  .  1.22 


Flqure  11.  Circumferential  Velocity  Profile!  for 
9  ■  0*.  90* .  130*.  270*  at  H  *  0.91 
a  -  2.0*.  X/0  -  S.5 


270°.  At  4  *  90°  cross  flow  velocity  caused  by  the  angle  of  attack  Is  In 
the  same  direction  as  the  wall  velocity.  At  4  *  270°,  the  outer  cross  flow 
velocity  opposes  the  wall  velocity.  The  circumferential  velocity  of  pro¬ 
files  at  ♦  *  0®  and  4  -  180°  are  equally  affected  by  the  surface  spin. 
Figure  11  shows  the  velocity  distribution  at  a  station  midway  on  the  boat- 
tall.  The  circumferential  velocity  distribution  at  4  »  90°  and  4  »  270° 
has  changed  significantly  from  that  shown  In  Figure  10.  On  the  boattall, 
the  decreasing  body  diameter  results  In  the  surface  velocity  decreasing  In 
magnitude.  However,  the  boundary  layer  thickness  In  this  region  Increases 
and  the  effect  of  surface  spin  Is  seen  to  persist  further  out. 
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CH  PLANE  AND  MAGNUS  COEFFICIENTS 


ntegration  of  the  pressure  and  viscous  forces  has  been  carried  out 

in  order  to  determine  the  aero¬ 
dynamic  coefficients.  The  sign 
convention  used  for  the  coeffi¬ 
cients  is  shown  in  Figure  12. 
The  results  in  Figures  13-17 
are  plotted  as  a  function  of 
longitudinal  position  and  thus 
show  how  the  force  develops 
over  the  length  of  the  projec¬ 
tile.  Figure  13  is  a  plot  of 
the  normal  force  coefficient 
and  shows  the  rapid  increase  in 
normal  force  which  occurs  on 
the  ogive  portion  of  the  pro¬ 
jectile.  The  cylinder  portion 
iroduce  no  significant  additional  normal  force;  however,  the  compu- 
ndicates  a  slight  increase  in  normal  force  here.  The  reversal  in 
iction  of  the  force  on  the  boattail  can  be  clearly  seen  as  the 
ited  normal  force  decreases  over  the  length  of  the  boattail.  The 
sntal  normal  force  coefficient,  indicated  by  the  circle,  shows  very 
•eement  with  the  computation. 


NORMAL  FORCE 


re  12.  Aerodynamic  Coefficient  Sign  Convention 


13.  Normal  Force  Coefficient  Along  the  Project) li 
Cowxjtatlon  and  Exoerlment,  H  .  0.91,  a  •  2. 


The  spin  rate  of  typical  artil¬ 
lery  shell  is  of  the  order  of 
300-500  rev/sec.  As  mentioned 
previously,  the  Magnus  effect 
results  from  a  spin  induced 
distortion  of  the  viscous  bound¬ 
ary  layer  which  occurs  for  artil¬ 
lery  shell  at  anale  of  attack. 
Previous  studies*’3  have  shown 
that  the  Magnus  effect  consists 
of  the  sum  of  the  boundary  layer 
displacement  effect  (asymmetric 
surface  pressure  distribution) 
plus  the  viscous  wall  shear 
«’  stress  contributions 


■gy)ymo  and  T$  “  w^y^y*o*  The  deve1°Pment  of  a11  three  components 

lagnus  force  are  shown  In  Figures  14,  15,  and  16  respectively,  as  a 
i  of  longitudinal  position.  Both  the  longitudinal  and  circumferen- 
jponents  (Figures  14  and  15)  are  seen  to  be  of  the  order  of  10" 5  and 
(spectlvely.  The  pressure  component  (Figure  16)  is  of  the  order 
rhe  dominant  component  of  the  transonic  Magnus  effect  is,  therefore 
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seen  to  be  the  boundary  layer 
displacement  effect,  Cyp. 

Additionally,  the  largest  portion 
of  the  total  Magnus  effect  is 
seen  to  develop  on  the  boattai? 
where  the  viscous  boundary  layer 
reaches  its  maximum  thickness. 
This  Is  the  same  qualitative 
behavior  reported  by  Sturek,  et 
al1’3  for  supersonic  flow.  The 
total  Magnus  force  (Cy  »  Cyu  + 

Cyw  +  Cyp)  Is  shown  In  Figure  17 

compared  to  the  experimental 
measurement.  Considering  the 
small  magnitude  of  the  Magnus 
force  and  the  agreement  achieved 
for  the  normal  force,  the  quanti¬ 
tative  agreement  between  the 
computation  and  experiment  is  regarded  as  very  good.  The  experimental 
Magnus  force  measurements  were  obtained  in  a  wind  tunnel  not  specifically 
designed  for  transonic  flow  and  are  considered  to  be  of  good  qualitative 
value  but  of  questionable  quantitative  value.  Additional  computations  at 
various  transonic  Mach  numbers  and  reliable  experimental  data  are  required 
before  a  full  assessment  of  the  computational  technique  can  be  made.  This 
first  result,  however,  for  predicting  the  Magnus  effect  at  transonic 
velocity  is  considered  very  encouraging. 


VII.  SUMMARY 

The  research  effort  presen¬ 
ted  in  this  paper  is  part  of  an 
overall  program  to  develop  a 
sophisticated  predictive  capabi¬ 
lity  for  projectile  aerodyna¬ 
mics.  The  pacing  requirement  for 
this  capability  is  the  determina¬ 
tion  of  the  Magnus  force  in  the 
transonic  flight  regime. 

An  implicit  finite  difference 
code,  which  solves  the  unsteady 
thin-layer  Navier-Stokes  equa- 


?tqur<  15.  LoKqltud'flai  MaP  '.h«ar  Contrioution  to  t*»  Naq*i'<«  tlOrtS,  h3S  been  applied  tO  3 

for»,  * •  o.«,  <> •  2.o*.  - .  333  r*v/**e  projectile  shape  at  a  ■  2.0° , 
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x/o 


M  ■  0.91.  The  solution  was 
marched  In  time  until  a  steady 
state  result  was  obtained. 
Computations  were  first  performed 
on  a  CDC  7600  using  a  finite  dif¬ 
ference  grid  of  21,600  points  and 
required  7.78  hours  of  computer 
time.  Increased  grid  resolution 
with  faster  computational  speed 
per  grid  point  was  obtained  by 
performing  the  computation  on  a 
Cray  IS  vector  computer. 

Tie  computations  have  been 
compared  to  experimental  surface 


Figur*  16.  circuaf«renti«i  Prt**ur«  contribution  to  tnt  *iqnus  pressures  and  aerodynamic  force 
Fore.  N  -  0.91.  .  .  i.v,  -  •  333  coefficients.  The  clrcumferen 


tlal  velocity  distribution,  presented  for  two  axial  locations,  showed  the 
significant  Interaction  between  the  cross  flow  velocity  resulting  from 
angle  of  attack  and  the  body  surface  velocity.  Experimental  velocity 
profile  data,  which  are  very  difficult  to  obtain  for  a  spinning  model  at 
transonic  speeds,  a»*e  required  to  fully  assess  the  computational  results. 

The  normal  and  Magnus  force  coef¬ 


ficients  have  been  shown  to  be  In 
good  quantitative  agreement  with 
experimental  data.  The  Individ¬ 
ual  components  of  the  Magnus 
force  have  additionally  been 
presented  and  Indicate  qualita¬ 
tively  good  results.  The  need 
for  additional  grid  resolution  or 
adaptive  grid  techniques11  have 
been  Identified  as  a  further 
requirement  to  achieve  more 
accurate  predictions.  Good 
quality  experimental,  transonic 
Magnus  data  is  also  required  for 
future  code  validation. 


The  present  results  Indicate  that  the  thin- layer  Navler-Stokes  computa¬ 
tional  technique,  in  conjunction  with  enhanced  computer  technology,  has  the 
potential  of  providing  the  capability  to  accurately  predict  the  aerodynamic 
behavior  of  spinning  shell  at  transonic  velocities,  Including  the  Magnus 
effect. 
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